Properties of quasi-one-dimensional molecules with Feshbach resonance interaction 
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Bound states and collisions of atoms with two-channel two-body interactions in harmonic waveg- 
uides are analyzed. The closed-channel contributions to two-atom bound states become dominant 
in the case of a weak resonance. At low energies and values of the non-resonant scattering length the 
problem can be approximated by a one-dimensional resonant model. Three-body problem becomes 
nonintegrable and the properties of triatomic molecules become different from those predicted by 
the integrable Lieb-Liniger-McGuire model. 

PACS numbers: 03.65.Ge, 03.65.Nk, 03.75.Be, 32.80.Pj 



(N 



O ■ 
>v 

43 ■ 



> 

m 
o 

o 

\o 
o 

o 

>>: 
43 : 

Oh- 



X 



Introduction 

Quasi-one-dimensional (ID) molecules, the existence 
of which was predicted in Rcf. 1], have been recently 
observed in an experiment Q with two-dimensional op- 
tical lattices. Except of lattices (see also Refs. 0, [3, p ) , 
tight cylindrical confinements, or atomic waveguides, 
have been realized in elongated atomic traps (see Refs. 
[E0 [S0 an d references therein), and atomic integrated 
optics devices (see Refs. and references therein). 

Unlike weakly bound molecules in free, three-dimensional 
(3D), space (see Refs. and references therein), ex- 

isting at positive elastic scattering length only, the quasi- 
1D ones survive at negative scattering length. 

The interaction of two atoms in an atomic waveguide 
can be considered as a ID zero-range interaction (see Ref. 
|14| ) whenever the collision or binding energies remain 
small compared to the transverse waveguide frequency 
u>x. Due to coupling to excited transverse states, the 
interaction strength demonstrates a resonant dependence 
on the ratio of the elastic scattering length a 3 r> to the 
transverse oscillator length 



aj_ 



(1) 



where m is the atomic mass. This confinement-induced 

resonance can been interpreted as a Feshbach resonance 
(see Refs. pj E3)> where the excited transverse states 
play the role of a closed channel. The interaction strength 
demonstrates also a dependence on the collision energy, 
which is a common property of a resonant scattering. 
Such an energy dependence appears in the scattering am- 
plitudes in Ref. djj] , as well as in the bound state energy 
in Ref. [Jj. However, the use of Feshbach resonance for 
tuning the elastic scattering length, as in Ref. 0, leads to 
additional energy dependence (see Ref. 0), which can 
be important for a weak resonance. A similar effect has 
been considered in Refs. |l9ll20ll2lL for a problem of 
two atoms under 3D harmonic confinement. 

A Feshbach resonance appears when the collision en- 
ergy of an atomic pair in the open channel lies in a vicin- 
ity of a bound (molecular) state in a closed channel (see 



Ref. H3)- As a result, the quasi- ID molecules are su- 
perpositions of the closed channel state and the ground 
and excited waveguide modes of the open channel. The 
closed-channel contribution becomes substantial for weak 
resonances. 

The theory of two-body two-channel problem under 
tight cylindrical harmonic confinement is summarized in 
Sec. [I] This problem can be approximated by a ID two- 
channel model. An improved relation between the ID 
and 3D scattering parameters, presented in Sec. [H] be- 
low, increases the range of applicability of the ID model 
compared to the relations in Ref. [18|. The composition 
of quasi-lD bound states is analyzed in Sec. IIIII Three- 
body ID molecules are considered in Sec. IIVI in a way 
which is similar to the analysis of scattering in Ref. [23| . 
A system of units with h = 1 is used below. 



I. FESHBACH RESONANCE IN HARMONIC 
WAVEGUIDES 

The properties of two-atom systems can be described 
by close-coupled equations for the wavefunction of the 
open channel ip a (r) and the amplitude for the system to 
be in the closed channel ip m , of the form (see Ref. |18|). 



--V 2 + U Q <5(r) + U conf (r) 
m 



(r) 



+V am S (r) ijj m 
Evb m = D 3DV j m + V a * m yj a (0). 



(2) 



Here E and r are, respectively, the energy and coordinate 
vector of the relative motion and all the energies are mea- 
sured from the open channel threshold. For a harmonic 
waveguide the confinement potential has the form 



T/ m 2 2 

Vcouf = 



(3) 



where p and z are the cylindrical components of the vec- 
tor r. The strength of the open channel potential V a , the 
coupling strength V am , and the bound state energy in the 
closed channel D^d can be expressed as (see Ref. p"8|) 



4tt 

V a = — a 3D 
m 



1 a3DP, 

7T 
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\V am \ 2 = —a 3D ^A (l - -a 3D p c 
m \ 7T 



(4) 



£>3_D = A* 



B-Bo-A + A 1 



-a3Z5Pf 



in terms of the phenomenological resonance strength 
A, the difference between the magnetic momenta of an 
atomic pair in the open and closed channels fi, the de- 
tuning of the external magnetic field B from its resonant 
value B , and the momentum cutoff p c . The final re- 
sults reached below are derived in the limit p c — ► oo. 
The (5-function potentials are applicable to two indistin- 
guishable bosons, as well as to bosons or fermions with 
different spins. 

Elimination of tjj m from Eqs. J3J) leads to a single equa- 
tion for ip a (r) . ft can be expanded in terms of the trans- 
verse Hamiltonian eigenfunctions |n0) with the zeroth 
angular momentum projection on the waveguide axis z 
as 



V>« (r) 



(2 7 r)- 1/2 E 



n=0 



dq^ n (q)e iflz \nQ). (5) 



The coefficients ip n (q) satisfy the set of coupled equations 
(see Ref. 0) 



m 



^ ' -ipn (q) 



oo 
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—^V cS (E)J2 / dq'^iq') 



where 



V eS (E) = V a + 



\V„, 



E- D 



(0) 
(7) 



is a non-renormalizcd energy-dependent interaction 
strength and 



Pn = y/m[E— (2n+ 1)cj±] 



(8) 



is the relative axial momentum for the channel corre- 
sponding to the transverse excited state |n0) with the 
excitation energy (2n + 1)lj±. 

The transition matrix for a two-atom collision in an 
atomic waveguide has been derived in Ref. 0] as 



T CO nf bo) 



ma i 



Ctefl (E) 



(9) 

Although the collision momentum p n depends on the 
channel, the transition matrix is independent of the ini- 
tial and final transverse states n, n' for all open channels 
(n,n' < (E/uj± — 1) /2) and is expressed by Eq. © in 
terms of po . It is a consequence of the use of zero-range 
potentials in Eqs. @. The energy-dependent length 



a e ft (E) = a 3D 



1 



E - n {B - B 



(10) 



replaces the elastic scattering length in Bethe-Peierls 
boundary condition. The Hurwitz zeta function is de- 
fined as (see Refs. Hill), 



C (y, a) — lim 



n c — 'oo 



n + a) 



1 



(n c + a) 



with — 2ir < arg(n + a) < 0. 



(11) 



II. RELATION TO THE ONE-DIMENSIONAL 
PROBLEM 

The confined two-body problem can be interpreted as 
a ID one described by the Schrodinger equation 



E c <Pa (z) = - 



1 d 2 tp 
m dz 2 



+ U eS (E c )5(z)<p Q (0) (12) 



with a zero-range interaction, where the interaction 
strength U c g depends on the collision energy E c = p 2 /m. 
The ID transition matrix corresponding to Eq. (II2II . 

Tid (po) = U, 



(13) 



coincides with Eq. © for 
1 



U c g(E c 



2u±a e g (E c 





E c \ 


im 




2co±) 


~ Wo. 



(14) 

The case of low collision energies E c <C 2ui i can be ana- 
lyzed using the expansion (see Refs. [l^Lllq') 



--C- C'a, ^\a\ 

a— >0 yja 



(15) 



1.3062, and 



where C = -( (§) « 1.4603, C = \Q (| 
£ (v) is the Riemann ^-function (see Ref. |24j). 

Following expressions attain a simpler form written in 
terms of dimensionless parameters (the scattering mo- 
mentum fc, the elastic scattering length a, the detuning 
b', and the resonance strength d), defined as 



b' = /x 



B-B -A 
2u\ 



2 ' 

1 

2' 



d = 



_ «3D 

2a±ujj_ 



(16) 



The par ameters fc, a, and d have been used previously in 
Ref. ■ The detuning 6', measured from the crossing 
point of the closed-channel bound state and the contin- 
uum threshold (D 3 q = uj±), is related to the detuning 
b of Ref. 0] , measured from the scattering length reso- 
nance, as 6' = b — d/a. 

Substitution of the expansion (I15|) into Eq. If 4|) leads 
to the expression of the interaction strength in the form 



t/cff (2u±k 2 



1/ 



ma ± C'ak 4 + (3k 2 - (I - Co) V 



d/a' 
(17) 
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where 



p = 1 - Ca - Cab'. 



(18) 



In a wide range of the parameters the interaction strength 
can be approximately expressed as 



U cS (E c ) = U a 



2|ff| : 



E c -D l 



D 



(19) 



This form corresponds to a two-channel ID problem (see 
Ref. 0]), described by the coupled equations 



Eipo (z) = - 



1 rfVo 
m dz 2 



U a po (0) + V2. 9 vH S(z) 



(20) 



E<p« m = D 1D tf m + V2gp Q (0) 



for the open- and closed-channel coefficients ipa (z) and 
ipi™ 1 , respectively. The non-resonant interaction strength 
U a , the channel coupling g, and the detuning Did will be 
further related to the 3D scattering parameters. Equa- 
tion (|T2*|) with U e s given by Eq. iJTSj l can be obtained by 
the elimination of the closed channel from Eqs. 1)20(1 . 

The ratio of the first term in the denominator of Eq. 
(|T7Jl to the other terms does not exceed the order of mag- 
nitude of ak 2 . It can be neglected whenever a < 1 and 
leading to an expression of the form of Eq. I(19|) 

with 



U, 



4o 



. |2 , d + C'{ab') 

\g\ = 4w -l 



M 2 



2w\ 



ma±(3 2 
1 - Ca) ab' + d 



a0 



(21) 



For the case of a relatively small detuning, C'\ab'\ <C 
|1 - Ca\, or 



\H (B - B - A - | < 



|0_lW_L 
C«3£> 



1 -C 



»3D 



the terms proportional to b' in the parameter (3 can 
be neglected. The parameters U a and Din&re ex- 
pressed then by Eqs. (43) and (45) in Ref. [II], while 
Eq. (44) therein will be valid whenever d ^> a 2 b' 2 , or 

p? (B - B - A - UJ±//J.) « a i a; ImA/qsdI. 

Unlike Eqs. (43)-(45) in Ref. the relations QTJ 

demonstrate a non-linear dependence of all three param- 
eters U a , g, and Did on the elastic scattering length, 
detuning, and resonance strength. These relations sub- 
stantially increase the applicability range of the ID ap- 
proximation, as is demonstrated by Fig. fusing the ex- 
ample of binding energy E u = u>± — E = 2u>±x 2 . The 
parameter x here is the ratio of a± to the bound state 
axial size. For the confined system it is determined as 
a solution of the transcendent equation (see Eq. (55) in 
Ref. [3) 



-u\ ( 1 x 2 j 0.. (2:!! 




FIG. 1: The binding energy Eb calculated as a function of the 
dimensionless resonance strength d [see Eq. llM l at the scat- 
tering length resonance B = Bq + cj±/fi, or 6' = —d/a, with 
the dimensionless non-resonant interaction strength a — 0.1. 
The solid, dashed, and dot-dashed lines present, respectively, 
solutions of the exact equation 123H . the ID approximation 
124H . and the open-channel model I25H . 



while for the related ID system it is evaluated as a solu- 
tion of the cubic equation 

f3x 3 + ax 2 + [(1 - Ca) b' + d/a] x + ab' = 0. (24) 

This equation determines poles of the ID T-matrix H13|) 
with the interaction strength l|19|) and the resonance pa- 
rameters given by Eq. 112111 . A similar equation has been 
considered in Ref. [25j ■ 

Substitution of the resonant scattering length 
&3D [1 — A/ (B — Bq — Lj_\_/fi)], which takes into account 
the resonance shift by u)±/fi due to confinement, into 
the equations of Ref. [l| , leads to the following equation 
for the parameter x 



ab' 



d Iff 1 2 



,x>0,. 



(25) 



This approximation, corresponding to a single-channel 
confined problem with an energy-independent interac- 
tion, is called here the 11 open- channel model' '. Figure 
H demonstrates that this model is applicable to strong 
resonances only. 

Various approximations for the binding energy are 
compared in Fig. |3 for a strong resonance. The re- 
sults demonstrate good agreement between solutions of 
the exact equation l)23|l . the ID approximation l)24[l. and 
the open-channel model 1)25(1 for B < Bq + A + 
when a weak bound state exists. However, the open- 
channel model predicts a non-physical singularity at B = 
B n + A + ui± j [i. Higher above the resonance the result of 
the open-channel model tends to the energy of the deep 
quasi-3D bound state. The latter state is not described 
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FIG. 2: The binding energy Eb calculated as a function of 
the external magnetic field for the 202 G resonance in K with 
u>± — 69 x 2-7T KHz. The solid line, pluses, and crosses present, 
respectively, solutions of the exact equation 12311 . the ID ap- 
proximation 1241 1. and the open-channel model 12511 . 
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FIG. 3: The binding energy Eb calculated as a function of 
the external magnetic field for the 543 G resonance in 6 Li 
with uj± = 200 x 27T KHz. The solid line, pluses, and crosses 
present, respectively, solutions of the exact equation il'i.'it . the 
ID approximation 1241 . and the open-channel model 125H . 



by the ID approximation. The exact equation l|23[l gives 
correctly both the deep and weak bound states. 

The ID approximation agrees to the exact equations 
also in the case of a weak resonance (see Fig- ■ However 
the open-channel model, neglecting energy dependence 
of the interaction strength, is applicable in this case only 
within a small interval of B close to Bq + A + cjx/m- 

The applicability range of the ID approximation is 
improved because the approximation (|15|) has a root at 
a = 0.294, close to exact value of 0.303 (cf. the value of 
0.468 provided by two-term approximation used in Ref. 



[l8|l. The accuracy of the binding energy demonstrates 
the applicability of the ID approximation to the off-shell 
T-matrix, which is important for applications to many- 
body problems (see also discussion in Ref. [26]). The 
only criteria of applicability are a < 1 and Eb <C uj±_ (or 
E c <C uj± for collisions). 

Equation 1(17)) can demonstrate an energy-dependence 
of the form of Eq. I|19|) for a weak resonance, whenever 
d <C ab', or 



A< \B - B - A - cjx/mI 



(26) 



and the term d/a in the denominator of Eq. I|17|) can 
be neglected. In this case the ID parameters can be 
estimated as 



U a « 0, \gf 



C 



-,D U 



-2lu 



1 - Ca 
C'a ' 



(27) 



They are independent of the Feshbach resonance pa- 
rameters. Therefore, in agreement with Ref. 0, a 
confinement-induced resonance can be interpreted as a 
two-state Feshbach resonance even for non-resonant 3D 
scattering. It is a consequence of the approximation 1)15(1 . 
It should be noted that the closed channel involves in this 
case a superposition of all transverse excitations. How- 
ever, the detuning Did then substantially exceeds the 
transverse frequency lu± and the energy-dependence is 
very weak in the quasi- ID regime, whenever the energy 
is less then u>±. 



III. BOUND STATE COMPOSITION 

The bound states of two atoms in atomic waveguides 
are superpositions of the closed and open channels of the 
Feshbach resonance. The open-channel component is a 
superposition of all transverse modes. The size of the 
closed-channel component is negligibly small compared 
to aj_, this component is not affected by the waveguide, 
and does not need an expansion in terms of the transverse 
modes. 

For a bound state with energy E < ujj_ all momenta 
p n , defined by Eq. (JSJ, are imaginary and the solution of 
Eq. 10 has the form 



$n {q) 



C a 



\Pn\ 



(28) 



The probability to find the molecule in the n-th trans- 
verse mode of the open channel can then be expressed 
as 



w = i!^ 



(29) 



Equations J2J) and JSJ allow to relate the closed-channel 
amplitude "0m to the open-channel wavefunction as 



0m = 



v* 

r am 



C a 



E - D 3D y/2- 



oo °9 

^2 / $ n ( q "> dq - 

n *f 



(30) 
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Although the sum here diverges, the renormalization pro- 
cedure of Ref. QJ3 leads to a finite probability to find the 
molecule in the closed channel 



1.0 



7r\C a \ 2 alnA 



2ma 3D [E - fi(B-B Q - A) 



(31) 



The final results can be expressed in terms of the dimen- 
sionless parameters 111 (il l 



2dx 3 



W a 2 (x 2 + b>) 2 ' 



We 

Wo 



x 3 ( 



2" 



1, (32) 



where x is the solution of Eq. igSJ, W e = J2 W n is the 

n=l 

total contribution of all excited transverse states, and the 
contributions are normalized as Wq + W e + W c — 1. An 
approximate expression, used in following calculations, 



l + 2x 2 



0.5 + 0.54884x + 2.5636a; 2 + 0.12172x 3 + 4x 4 



gives a relative error of less then 3 x 10~ 4 . 

Bound states of the related ID system are superposi- 
tions of open and closed channels, where the closed chan- 
nel effectively incorporates contributions of the excited 
transverse states and the closed channel of the confined 
system. The two contributions can be respectively ex- 
pressed as 



Wl D = \cpf 



l^o (z) \ 2 dz, (33) 



where the bound-state solution of Eqs. I|2U|) has the form 
V2g<p (0) 



E-D 



— , ipo (z) = ip (0) exp (~\JmE b \z\ 
id v 



The ratio of the contributions can be expressed as 
Wl D . , d + C'a 2 b 12 



= 2x i 



Wl D a 2 (x 2 + b>) 2 ' 

A direct evaluation demonstrates that 

W ° + W < ~ ^ where C (I A - 



(34) 



(35) 



(36) 



The contributions to the bound states are presented in 
Fig. El for two resonances: the strong one in K with ad « 
1.6 and the weak one in 6 Li with ad ps 4.1 x 10~ 4 . The 
results demonstrate that W} D &W e + W c and W^ D m 
Wq with a good accuracy over a wide range of parameter 
values. The closed channel and excited states of the open 
channel yield the dominant contribution far below the 
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FIG. 4: (color online) The contributions of the ground (dot- 
dashed line) and excited (dashed line) states of the open chan- 
nel and of the closed channel (solid line) to the weak bound 
state calculated as a function of the external magnetic field. 
The dot-dot-dot-dashed and dotted lines present, respectively, 
contributions of the open and closed channels in the ID model. 
The parts (a) and (b) are related, respectively, to the 202 G 
resonance in K with uj± = 69 x 2n KHz and to the 543 G 
resonance in 6 Li with u>j_ — 200 x 2ir KHz. The dot-dashed 
and dot-dot-dot-dashed lines almost coincide in the part (b). 



bound state threshold at B = Bo + A + (J±/[i, while the 
ground state of the open channel becomes dominant near 
the threshold. For the strong resonance the contribution 
of the closed channel in the confined system is always 
small, and the bound state consists mostly of the ground 
and excited states of the open channel. However, for 
the weak resonance the major contribution is yielded by 
the closed channel and the ground state of the excited 
channel. The closed channel and excited states of the 
open channel yield similar contributions far below the 



bound state threshold for ad — ma 2 B fiA/4: w 0.1. The 
resonances with higher or lower values of the product 
ad can be called, respectively, as open-channel or closed- 
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channel dominated resonances, as in a case of free space 
(see Refs. 003). 



IV. THREE-ATOM BOUND STATES 

The previous results demonstrate that two atoms with 
a resonant interaction in an atomic waveguide can be 
within a good accuracy considered as ID particles with 
a resonant interaction. Consider now a three-body ID 
problem for bosonic atoms. (The same approach has 
been used in Ref. [23| for the analysis of three-body scat- 
tering.) A state vector can be represented in the form 



I*,) = 



dzdz m ipf ) (z, z m ) i>\ (z) & m (z m ) 




V6 



d 3 Z^ ] ( ZU Z2,Z 3 )& a ( Zl )¥ a (z 2 )*t (z 3 ) 



10*37) 



as a superposition of the three-atom channel, described 
by the coefficient (p (zi, z 2 , z 3 ), and the atom-molecule 
channel (involving the closed-channel molecules), de- 
scribed by the coefficient ip^' [z,z m ). Here ^\ (z) and 
^ m (z m ) are the creation operators for the atom and 
closed-channel molecule, respectively, and z, z m are 
their coordinates. Substitution of Eq. lpF7|) into the 
Schrodinger equation with the Hamiltonian (25) of Ref. 
[l8| leads to the following coupled equations 



Eip { 3) (Z!,Z 2 ,Z 3 ) = 



1 3 

— y 

~>m 



2m ' — ' Cz„ 

j=i j 



. . ' i U a [5 (zi - z 2 ) 
az z 



FIG. 5: The binding energy Eb calculated as a function 
of the dimensionless detuning (see Eq. 1421 '). The solid, 
dashed, and dot-dashed lines present the energies of three- 
body bound states for the dimensionless non-resonant inter- 
action strengths u = 1,0, and -1, respectively. The corre- 
sponding energies of the two-body bound states are presented 
by pluses, circles, and crosses, respectively. 



x<4 3) (z!,z 2 ,Z3)(39) 
obeys the single-channel Schrodinger equation 

1 3 

E <fo (Qi,Q2, 93) = -z— ^2 ? f ™ (91^2, 93) 



+6 (z 2 - z 3 ) + 5 (zi - z 3 ) 



(3) I \ 
</?0 (Zl,22,Z 3 ) 



+<y5i 3) (z 2 , zi) 5 (zi - z 3 ) + <p { i' (z 3 , zi) 5 [zi - z 2 

id 2 id 2 



+ \l ^9* [f i 0*1, 22) 6 (z 2 - z 3 ) 



3=1 



Etp^ (z, z m ) = 



x / d 3 q'S {q' 3 - Qj ) S (Q - Q') Cpf q' 3 ) . (40) 



2m dz 2 Am dz 2 



D 



(f^> (Z, Zm) 



The indistinguishability of the bosonic atoms leads to 
the symmetry of the wavefunction ip (zi,z 2 ,z 3 ) over 
permutation of the atomic coordinates. 

Equations (|38l) . as well as the Hamiltonian, do not 
contain terms describing collisions between the closed- 
channel molecule and the third atom. This assumption 
is justified since the atoms in the closed and open chan- 
nels have different hyperfine states, ft allows a simple 
elimination of the atom-molecule channel. The analysis 
of the remaining three-atom channel attains a simpler 
form in the momentum representation. The correspond- 
ing three-atom wavefunction 



Here qj are the atomic momenta and Q = q% + q 2 + q 3 is 
the total momentum. The interaction strength U Q g hap- 
pens to be the same function [see Eq. ifHEjl ] as in two-body 
+V6g<p Q (z, z m , z m ) problem. The conventional Faddeev reduction technique 
(see Ref. [27j ) leads for three-body bound states to the 
homogeneous equation 



X(q) = ^ J dq' 



mE — q' 2 — q'q — q 2 



l fio ) (91,93,93) 



(27T) 



-3/2 



d gexp 




<T U > I i K m\E\ + -<f* \X(q'), (41) 



where E < is the three-body bound state energy in 
the center-of-mass system and the ID transition matrix 
is given by Eq. (|T3|) . 

Properties of ID systems depend on two dimensionless 
parameters: the non-resonant interaction strength and 
the detuning, respectively, 

u = m^\g\- 2 / 3 U a , b 1D = m-^\g\-^D w . (42) 
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Figure presents the scaled binding energy e 3 = 
_ m -i/3|g|-4/3£; ca i cu i a ted by a numerical solution of 
Eq. 141(1 . A related problem with non-resonant interac- 
tions only, the Lieb-Liniger-McGuire model (2^, [2£| , has 
an exact solution. The binding energies for two- and 
three-body bound states in that model are expressed as 
e 2 = —u 2 /2 and e 3 = —2m 2 , respectively, and therefore 
£3 = 4e 2 . In the present resonant case £3 ~ 4e 2 only 
at large positive detunings. For large negative detun- 
ings the two-body bound state contains mostly the ID 
closed-channel contribution, and the three-body bound 
state has a form of the two-body state with a third atom 
weakly bounded to it. 

Conclusions 

Two atoms with a Feshbach resonant interaction in 
an atomic waveguide can be approximated by a ID reso- 
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